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Abstract 

RNA-dependent RNA polymerases (RDRs) act as key components of the small RNA biogenesis pathways and play significant 
roles in post-transcriptional gene silencing (PTGS) and antiviral defense. However, there is no information about the RDR 
gene family in Salvia miltiorrhiza, an emerging model medicinal plant with great economic value. Through genome-wide 
predication and subsequent molecular cloning, five full-length 5. miltiorrliiza RDR genes, termed SmRDRl-SmRDRS, were 
identified. The length of SmRDR cDNAs varies between 3,262 {SmRDRS) and 4,130 bp {SmRDRS). The intron number of 
SmRDR genes varies from 3 {SmRDRl SmRDRS and SmRDR4) to 17 (SmRDRS). All of the deduced SmRDR protein sequences 
contain the conserved RdRp domain. Moreover, SmRDR2 and SmRDR4 have an additional RRM domain. Based on the 
phylogenetic tree constructed with sixteen RDRs from Arabidopsis, rice and 5. miltiorrliiza, plant RDRs may be divided into 
four groups (RDR1-RDR4). The RDR1 group contains an AtRDR and an OsRDR, while includes two SmRDRs. On the contrary, 
the RDR3 group contains three AtRDRs and two OsRDRs, but has only one SmRDR. SmRDRs were differentially expressed in 
flowers, leaves, stems and roots of 5. miltiorrhiza and responsive to methyl jasmonate treatment and cucumber mosaic virus 
infection. The results suggest the involvement of RDRs in 5. miltiorrhiza development and response to abiotic and biotic 
stresses. It provides a foundation for further studying the regulation and biological functions of SmRDRs and the biogenesis 
pathways of small RNAs in 5. miltiorrhiza. 
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Introduction 

RNA-dependent RNA polymerases (RDRs), which catalyze the 
conversion of single-stranded RNAs (ssRNAs) into double- 
stranded ones, play vital roles in the production of various small 
interfering RNA (siRNA) species in plants through collaboration 
with other proteins, such as Dicer-likes (DCLs) capable of cleaving 
double-stranded RNAs (dsRNAs) into 21-24 nt duplexes [1,2]. 
Based on the origins and biogenesis pathways, siRNAs generated 
from dsRNAs can be classified into several groups, such as trans- 
acting small interfering RNAs (ta-siRNAs), heterochromatic 
siRNAs (hsiRNAs) and natural antisense transcript-derived 
siRNAs (nat-siRNAs) [3-5] . These siRNAs may incorporated into 
the Argonaute (AGO)-containing RNA-induced silencing com- 
plexes (RISCs) to silence a variety of gene transcripts, repetitive 
sequences, sense transgenes, viruses and mobile elements through 
RNA cleavage, translational inhibition, DNA methylation and 
heterochromatin formation [6]. 

RDR proteins are characterized by the conserved RNA- 
dependent RNA polymerase catalytic domain (RdRp) and are 
among the first components identified for plant small RNA 
biogenesis pathways [7-10]. They are present in fungi, viruses. 



plants and nematodes, but have not been found in insects and 
vertebrates [11,12]. The activity of RDR was detected in Chinese 
cabbage more than forty years ago [8], whereas the cDNA of RDR 
was first isolated from tomato in 1998 [13]. So far, RDR genes 
have been identified in various plant species, such as Arabidopsis 
thaliana, rice, maize, tomato and tobacco [14-17]. Similar to DCLs 
and AGOs, RDR genes exist as a family in plants and the number of 
RDR genes may be not the same in different species. For instance, 
there are six RDR members in Arabidopsis and Solanum lycopersicum 
[14], five in rice [15] and maize [16], and at least three in Nicotiana 
tobaccum and N. attenuate [17]. Each small RNA biogenesis pathway 
may involve different member of the RDR family. 

The six members included in the A. thaliana AtRDR gene family 
were termed AtRDR l-AtRDR6, respectively [14]. AtRDRl, 
AtRDR2 and AtRDRG play distinct and overlapping functions in 
various aspects, such as viral resistance, chromatin silencing and 
post-transcriptional gene silencing (PTGS) [19,20]. AtRDRl and its 
ortholog in tobacco, MRDRl, are induced by salicylic acid (SA) 
treatment and virus infection and involved in plant susceptibility to 
tobacco mosaic tobamovirus (TMV) and tobacco rattle virus 
(TRV) [2 1] . The underlying functional mechanism of AtRDRl 
appears to produce and amplify exogenous and virus-derived 
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Table 1. Sequence features and intron number of SmRDRs. 







Gene name 


Accession number 


cDNA (bp) 


ORF (bp) 


5'UTR (bp) 


3'UTR (bp) 


Protein (aa) 


Mw (kDa) 


p/ 


Intron no. 




SmRDRI 


KF872203 


3606 


3357 


57 


192 


1118 


127.93 


8.22 


3 


SmRDR2 


KF872204 


3609 


3306 


148 


155 


1101 


125.55 


7.58 


4 


SmRDR3 


KF872205 


4130 


3588 


293 


249 


1195 


136.33 


7.67 


3 


SmRDRA 


KF872206 


3537 


2919 


35 


583 


972 


110.04 


6.92 


3 


SmRDRS 


KF872207 


3262 


2697 


65 


500 


898 


102.53 


7.07 


17 





doi:1 0.1 371/joumal.pone.00951 1 7.t001 



siRNAs (vsiRNAs) in infected plants [19-24]. In addition to 
antiviral responses, RDRl plays a significant role in plant 
resistance to herbivore attack [18]. AtRDR2 is involved in the 
production of the most abundant endogenous hsiRNAs that are 
mostly 24 nt in length and are associated with heterochromatic 
and repetitive regions, such as the pericentromeric regions and 
telomeres [25]. RDR2-dependent RNA-directed DNA methyla- 
tion (RdDM) is responsible for siRNA-mediated DNA methylation 
and histone modifications at Arabidopsis telomeres, and is required 
for the maintenance of telomeric heterochromatin [25]. In 
addition, AtRDR2 is involved in the development of the female 
gametophyte [26]. AtRDR6, acts in various gene silencing 
pathways and plays important roles in the biogenesis of ta-siRNA 
and nat-siRNA. It is also well-known for the amplification of 
improper terminated and unpolyadenylated RNAs generated from 
transgenes or inverted repeats to trigger degradation of comple- 
mentary RNA species [27]. The functions of AtRDR3~AtRDR5 are 
currently unknown [7]. 

Although RDR genes have been isolated from various plant 
species, to our best knowledge, there is no report on RDRs in 
medicinal plants. Salvia miltiorrhiza, well-known as danshen in 
Chinese, is an economically significant medicinal plant and is an 
emerging model medicinal plant for Traditional Chinese Medicine 
(TCM) studies [28]. It has been used for treating various human 
diseases, such as dysmenorrhoea, amenorrhoea and cardiovascular 
disease, for thousands of years [28-30]. The ^S*. miltiorrhiza genome 
has been preliminarily decoded (Chen et al, unpublished data). 
We have previously analyzed the S. miltiorrhiza AGO gene families 
involved in small RNA production and action [31]. In order to 
characterize the RDR genes in ^S*. miltiorrhiza^ we performed a 
genome-wide search of RDRs against the working draft of the S. 
miltiorrhiza genome, followed by molecular cloning and molecular 
analysis. The results will be a basis for understanding the gene 
silencing pathways in S. miltiorrhiza. 

Materials and Methods 

Plant materials and stress treatment 

S. miltiorrhiza Bunge (line 993) were grown in a field nursery. 
Flowers, leaves, stems and roots were collected from two-year-old 
plants and stored in liquid nitrogen until use. Plantlets cultivated in 
vitro were grown at 25 °C with a photoperiod of 16 h light and 8 h 
dark for six weeks as described previously [32]. MeJA treatment 
were carried out following the procedures reported previously 
[32,33]. Plantlets were treated for 12, 24, 36 and 48 h and then 
sampled. Sterile water-treated plantlets were used as controls. For 
cucumber mosaic virus (CMV) infection, the silicon carbide 
powder friction method was used. Briefly, leaves of six-week-old 
plantlets cultivated in vitro were dusted with silicon carbide powder 



and then inoculated with CMV subgroup I for 12, 24, 48 and 
72 h. Leaves inoculated with phosphate buffered saline (PBS) were 
used as controls. All tissues collected were stored in liquid nitrogen 
until use. Three independent biological replicates were performed 
for each experiment. 

Identification of SmRDR genes 

Arabidopsis and rice RDR protein sequences were downloaded 
from GenBank (http://www.ncbi.nlm.nih.gov/protein). It in- 
cludes AtRDRl (AEE29226.1), AtRDR2 (AEE82976.1), AtRDR3 
(082 190.2), AtRDR4 (082 189.2), AtRDR5 (082 188.2), AtRDR6 
(AEE78550.1), OsRDRl (QpDXS3.2), OsRDR2 (Q,7XM31.1), 
OsRDR3 (Q5QMN5.2), OsRDR4 (Q5QMN4.2), and OsRDR6 
(Q8LHH9.1). S. miltiorrhiza SmRDR genes were predicted by 
BLAST analysis of Arabidopsis and rice RDRs against the working 
draft of the S. miltiorrhiza genome (Chen et al., unpublished data) 
using tBLASTn [34,35]. The retrieved genomic DNA sequences 
were used for gene model prediction on the GENSCAN web 
server (http://genes.mit.edu/GENSCAN.html). Gene models 
were manually corrected according to the alignment between 
SmRDRs and other plant RDRs obtained from BLAST analysis 
of predicted SmRDRs against the non-redundant protein 
sequence (nr) database using the BLASTx algorithm (http:// 
www.ncbi.nlm.nih.gov/BLAST). 

Molecular cloning of SmRDR cDNAs 

The 5 '-RACE and 3 '-RACE were carried out as described 
previously [31]. Briefly, total RNA extracted from the root of S. 
miltiorrhiza was purified using the oligotex mRNA mini kit 
(Invitrogen). 5' and 3' RACE was performed on mRNA using 
the GeneRacer kit (Invitrogen). Gene specific nesting and nested 
primers were designed and synthesized (Tables SI and S2). PGR 
products were purified, cloned and sequenced. Based on the 
obtained 5' and 3' cDNA sequence, gene-specific forward and 
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Figure 1. Gene structures of 5. miltiorrhiza SmRDRs, Filled boxes 
represent exons with coding regions in green and 5'- and 3'-UTRs in 
blue. The connecting lines represent introns. 
doi:10.1371/journal.pone.0095117.g001 
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Figure 2. Conserved domains in SmRDR proteins. Boxes with 
points represent the RRM domain. Grey boxes represent the RdRp 
domain. The first and the last amino acids of SmRDR proteins and the 
location of RRM and RdRp domains are numbered. 
doi:10.1371/journal.pone.0095117.g002 

reverse primers were designed and synthesized (Table S3). Full- 
length SmRDR cDNAs were PCR-amplified, cloned and se- 
quenced as described [31]. 

Bioinformatic analysis and phylogenetic tree 
construction 

Bioinformatic analysis of SmRDR sequence features, such as 
intron/exon structures, molecular weight (MW), theoretical 
isoelectric point (pi) and conserved domain, were performed as 
described previously [31,36]. The conserved motifs of SmRDRs 
were analyzed using Multiple Expectation Maximization for Motif 
Elicitation (MEME) version 4.9.1 [37] with the following 
parameters. Optimum motif width was set to>6 and ^50. The 
maximum number of motifs was designated to identify 20 motifs 
following the previous reported studies [17]. The conserved 
residues were analyzed by alignment of amino acid sequences 
using T-coffee [38]. Phylogenetic tree was constructed for protein 
sequences of sixteen RDRs from S. miltiorrhiza, Arabidopsis and rice 
using MEGA version 4.0 by the neighbor-joining method with 
1,000 bootstrap replicates [39,40]. 



and CMV was analyzed using the qRT-PCR method as described 
previously [31]. Briefly, Gene-specific forward and reverse primers 
were designed and synthesized (Table S4). About 10 ng cDNA 
reversely transcribed from total RNA was used as a template in a 
20 |il volume. SmUBQlO was used as a reference [31]. qPCR was 
carried out in triplicates for each biological sample using the BIO- 
RAD CFX system (Bio-Rad). Three fully independent biological 
replicates were performed. The specificity of amplification was 
assessed by dissociation curve analysis. Gene expression levels were 
determined using the 2 ^^^^ method, where Cq represents the 
threshold cycle [41]. Relative amount of transcripts was calculated 
and normalized as described previously [42]. Average Cq were log 
transformed, mean centered and autoscaled [42]. Standard 
deviations of mean value from three biological replicates were 
calculated as described previously [42]. 

Results 

Identification and molecular cloning of five 5. miltiorrhiza 
RDR genes 

BLAST analysis of Arabidopsis and rice RDRs against the 
working draft of the S. miltiorrhiza genome (Chen et al., 
unpublished data) using tBLASTn [34,35] showed the existence 
of five SmRDR gene loci in the S. miltiorrhiza genome. Genomic 
DNA sequences were retrieved and predicted for gene models on 
the GENSCAN web server (http://genes.mit.edu/GENSCAN. 
html). The five gene models computationally predicted were 
BLAST-analyzed against the non-redundant protein sequence (nr) 
database (http://www.ncbi.nlm.nih.gov/BLAST) using BLASTx 
with default parameters and then manually corrected according to 
the alignment between SmRDRs and other plant RDRs. To 
further experimentally validate the predicted cDNA sequences of 
SmRDRs, molecular cloning of full-length SmRDR cDNA was 
carried out using RNA ligase-mediated rapid amplification of 5 ' 
(5' RACE) and 3' (3' RACE) cDNA ends and subsequent PGR 
amplification of coding regions. The deduced amino acid 
sequences of all five SmRDRs share high sequence identity with 
known plant RDRs and contain the conserved RdRp domain, 
suggesting they are authentic SmRDRs. The identifed SmRDR 
genes are termed SmRDR l-SmRDR5, respectively. The cloned 
cDNAs have been submitted to GenBank under the accession 
numbers KF872203-KF872207. 



Quantitative real-time reverse transcription-PCR (qRT- 
PCR) 

Expression of SmRDRs in roots, stems, leaves and flowers of 2- 
year-old ^S*. miltiorrhiza plants and in plantlets treated with MeJA 



Gene structure and conserved domain analyses 

Sequence feature analysis of SmRDRs suggests that the length of 
open reading frames (ORFs) of SmRDRs varies from 2,697 
{SmRDRS) to 3,588 bp {SmRDRS) (Table 1). The length of 5' and 
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Figure 3. Conserved motifs of SmRDRs proteins identified with the MEME search tool. Motifs are represented by boxes. The numbers (1- 
20) and different colors in boxes represent motif 1 -motif 20, respectively. Box size indicates the length of motifs. Broken lines indicate locations of the 
RdRp domain. 

doi:10.1371/journal.pone.0095117.g003 
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Figure 4. Alignment of partial sequences of RDR proteins in 5. miltiorrhiza and Arabidopsis using T-coffee. The DLDGD/DFDGD signature 
is boxed. Different colors represent the quality of alignment with red for the highest quality, yellow for average, and green for the worst. 
doi:10.1371/journal.pone.0095117.g004 



3' UTRs varies between 35 and 293 bp and between 155 and 
583 bp, respectively. The size of deduced SmRDR proteins varies 
between 898 and 1195 amino acids. The molecular weight (Mw) 
varies from 102.53 to 136.33 kDa, and the theorectical p/ is 
between 6.92 and 8.22 (Table 1). The SmRDR5 locus, which 
produces the shortest SmRDR, has seventeen introns (Fig. 1, 
Table 1). SmRDR2 contains four introns including one is located 
within the 5' untranslated region (UTR) (Fig. 1, Table 1). The rest 
three, including SmRDRl, SmRDRS and SmRDR4, have 3 introns 
(Fig. 1, Table 1). All of the introns of SmRDRl and SmRDR4 are 
located in the coding regions; however, of the 3 introns of 
SmRDRS, only one is located in the coding region. The other two 
are located within the 5' UTR. These 5' UTR-located introns 
might enhance gene transcription and RNA stability [43-45]. 
Additionally, the size of introns varies significantly among SmRDRs 
(Fig. 1). It suggests the diversity of SmRDRs in sequence features 
and gene structures. 

The search for conserved domains in SmRDR proteins against 
the NGBI Conserved Domain Database showed that all of the five 
SmRDRs contained the conserved RdRp domain (Fig. 2). It is 



consistent with the results from other plant RDRs [15-17]. 
SmRDR2 and SmRDR4 have an additional RNA recognition 
motif (RRM) in the region close to the N-terminus (Fig. 2). RRM 
in SmRDR2 starts from the 4* amino acid and ends at the 66*, 
while SmRDR4 RRM starts from the 1 4* amino acid and ends at 
the 9V\ The actual function of additional RRM in RDR proteins 
remains to be elucidated [46]. 

Using the MEME motif search tool, we analyzed conserved 
motifs of SmRDRs (Fig. 3). The results revealed four motifs (1, 7, 
10 and 12) conserved in all of the five SmRDRs. It suggests the 
conservation of SmRDRs. On the other hand, various less 
conserved motifs were found. For instance, motifs 11, 15, 16, 17 
and 19 are specific to SmRDRl and SmRDR2. Motifs 13 and 20 
are specific to SmRDRl, SmRDR2 and SmRDR3. Additionally, 
among the five SmRDRs, SmRDRl and SmRDR2 have 20 
motifs, whereas SmRDR5 only contain 4. These less conserved 
motifs could be associated with gene-specific functions. Sequence 
alignment of RDR proteins from S. miltiorrhiza and Arabidopsis using 
T-coffee [38] showed that SmRDRs contained the DLDGD/ 
DFDGD signature, the partial sequence of motif 1 (Figs. 3 and 4). 
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Figure 5. Phylogenetic relationships of sixteen RDRs from 5. miltiorrhiza, Arabidopsis and rice. The relationships were analyzed for 
deduced full-length amino acid sequences using MEGA 4.0 by the neighbor-joining (NJ) method with 1000 bootstrap replicates. Bootstrap values are 
shown near the nodes. Four groups of RDRs, termed RDR1 , RDR2, RDR3 and RDR4, respectively, are indicated. The number of introns in open reading 
frame (ORF) of the corresponding RDR gene is shown in parentheses. 
doi:10.1371/journal.pone.0095117.g005 
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Figure 6. Expression of SmRDRs in flowers (Fl), leaves (Le), 
stems (St) and roots (Rt) of 5. miltiorrhiza. The expression patterns 
were analyzed using the quantitative RT-PCR method. PCR was carried 
out in triplicates for each biological sample. Three independent 
biological replicates were performed. SmUBQIO was used as a 
reference. Fold changes of SmRDR expression are shown. The levels in 
roots were arbitrarily set to 1 and the levels in other tissues were given 
relative to this. Error bars represent the standard deviations of the mean 
value of three biological replicates. 
doi:10.1371/journal.pone.0095117.g006 

DLDGD/DFDGD has been previously found in various other 
plant RDRs and seems to be part of the nucleotidyl transferase 
active site of RDR proteins [47] . It further confirms the role of 
identified SmRDRs in the conversion of ssRNAs into dsRNAs. 

Phylogenetic analysis of RDR proteins in S. miltiorrhiza, 
Arabidopsis and rice 

Previous study revealed that RDRs in eukaryotic organisms 
might be divided into three clades, RDRcx, RDR(3, and RDRy 
[47]. RDRocproteins exist in all three kingdoms, whereas the 
proteins included in the RDRP clades are present only in animals 
and fungi and RDRyproteins are found only in plants and fungi. 
Among the six Arabidopsis AtRDRs, AtRDRl, AtRDR2 and 
AtRDR6 belong to the RDRa clade, while AtRDR3, AtRDR4 
and AtRDR5 are included in the RDRyclade [47] . To determine 
the evolutionary relationship among RDRs from S. miltiorrhiza, 
Arabidopsis and rice, an unrooted neighbor-joining tree was 
constructed for the full-length protein sequences of five SmRDRs, 



six AtRDRs and five OsRDRs. The results showed that the sixteen 
RDRs might be divided into four groups, termed RDRl, RDR2, 
RDRS and RDR4, respectively (Fig. 5). It is consistent with 
previous results for RDRs from Arabidopsis, rice and ^ea mays [15- 
17]. The RDRl, RDR2 and RDR4 groups contain SmRDRI, 
SmRDR2, SmRDRS and SmRDR4 from S. miltiorrhiz^, AtRDRl, 
AtRDR2 and AtRDR6 from Arabidopsis, and OsRDRl, OsRDR2 
and OsRDRG from rice. All of them belong to the RDRocclade 
and are characterized by the DLDGD signature (Fig. 4). 
SmRDRI and SmRDR2 and the antiviral defense-associated 
AtRDRl and OsRDRl [15,16] are included in the RDRl group. 
Members of the RDR2 group includes SmRDR4, Arabidopsis 
AtRDR2 and rice OsRDR2, of which AtRDR2 and OsRDR2 are 
involved in the production of the most abundant endogenous 
hsiRNAs and are responsible for siRNA-mediated DNA methyl- 
ation and histone modifications at telomeres [25]. The RDR4 
group contains SmRDRS, AtRDR6 and OsRDRG. AtRDRG is 
associated with amplifying improper terminated and unpolyade- 
nylated RNAs generated from transgenes or inverted repeats [27]. 
The RDRS group includes AtRDRS, AtRDR4 and AtRDR5 from 
Arabidopsis, OsRDRS and OsRDR4 from rice and SmRDR5 from 
S. miltiorrhiza. Members of the RDRS group belong to the 
RDRyclade and are characterized by the DFDGD signature 
(Fig. 4). Although it is the biggest group, the function of RDRs in 
the RDRS group is currently unknown. 

Tissue-specific expression of SmRDR genes 

To preliminarily elucidate the function of SmRDR genes, we 
analyzed the expression patterns of five identified SmRDRs in 
flowers, leaves, stems and roots of 2-year-old and field nursery- 
grown S. miltiorrhiza using the quantitative RT-PCR technology. 
The transcripts of all five SmRDRs could be detected in the tissues 
analyzed (Fig. 6), which is consistent with the vital roles of RDRs in 
plants. All of them showed the highest expression in roots and less 
in flowers and leaves (Fig. 6). Further sequencing and analyzing the 
sRNAome in ^S". miltiorrhiza may help to elucidate the underlying 
mechanisms. 

The response of SmRDRs to MeJA treatment and CMV 
infection 

RDRs catalyze the conversion of single-stranded RNAs into 
double-stranded ones and are core components for the production 
of siRNAs involved in plant development and response to abiotic 
and biotic stresses. In order to investigate the expression pattern of 
SmRDR genes under abiotic treatments, the expression level of 
SmRDRs in leaves of plantlets treated with MeJA was analyzed 
using the quantitative RT-PCR method. Plantlets treated with 
sterile water were used as controls. The results showed that all five 
SmRDRs were suppressed under MeJA treatment (Fig. 7). MeJA 
usually up-regulates the expression of genes associated with the 
biosynthesis of secondary metabolites, which play significant roles 
in plant response to stress [S2,SS]. Down-regulation of SmRDR 
gene expression could be helpful to increase the production of 
stress-related secondary metabolites. 

To examine the response of SmRDRs in biotic stress, S. 
miltiorrhiza plantlets were inoculated with CMV using the silicon 
carbide powder friction method. The level of SmRDR transcripts in 
leaves treated for 12, 24, 48 and 72 h was analyzed using the 
quantitative RT-PCR method. Leaves inoculated with phosphate 
buffered saline (PBS) were used as controls. As shown in Fig. 8, 
SmRDRI, SmRDR2 and SmRDRS were up-regulated at different 
time points. The level of SmRDRI showed a 1 .9-fold-increase after 
CMV inoculation for 48 h (Fig. 8A). SmRDR2 showed 2.S-, S.4- 
and 4.2-fold-increase after being treated with CMV for 12, 48 and 
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Figure 7. SmRDRs responsive to MeJA treatment. The expression patterns were analyzed using the quantitative RT-PCR method. PCR was 
carried out in triplicates for each biological sample. Three independent biological replicates were performed. SmilBQW was used as a reference. Fold 
changes of SmRDRs in leaves of 5. miltiorrhiza plantlets treated with MeJA for 12, 24, 36 and 48 h are shown. The level of transcripts in leaves treated 
with sterile water (CK) was arbitrarily set to 1 and the level in tissues treated with MeJA was given relative to this. Error bars represent standard 
deviations of mean value from three biological replicates. 
doi:1 0.1 371 /journal. pone.0095117.g007 



72 h, respectively (Fig. 8B). SmRDRS was accumulated to 3.7- and 
2.3-folds of controls after CMV infection for 12 and 72 h (Fig. 8C). 
However, no significant change was observed for the level of 
SmRDR4 and SmRDR5 after CMV infection (Figs. 8D and BE). 
The results indicate that SmRDRI, SmRDR2 and SmRDRS may be 
involved in antiviral defense in S. miltiorrhiza. 

Discussion 

S. miltiorrhiza is an economically significant medicinal plant 
species belonging to the largest genus, Salvia, in the mint family. It 
is native to China and Japan and is widely distributed in China. S. 
miltiorrhiza has close phylogenetic relationships with other Asian 
and Mediterranean Salvia species, such as S. roborowskii and S. 
glutinosa [48]. The root of ^S". miltiorrhiza has been widely used in 
TCMs for hundred of years to treat dysmenorrhoea, amenorrhoea 
and cardiovascular diseases [28-30]. The main bioactive compo- 
nents in S. miltiorrhiza are lipophilic diterpenoid tanshinones and 
hydrophilic phenolic acids. Genes involved in the biosynthesis of 
these components have been intensely studied recently [32,33,49- 
5 1] . Because of its relatively small genome size, short life cycle, 
undemanding growth requirements and significant medicinal 
value, S. miltiorrhiza is being developed to become a model 
medicinal plant for TCM studies [31-33]. Elucidation of small 
RNA biogenesis pathways in S. miltiorrhiza appears to be urgent, 
given the significant regulatory roles of small RNAs in plant 
development and growth. Results from Arabidopsis suggest that the 



core components of small RNA pathways include at least three 
gene families, DCL, AGO and RDR [2]. The AGO gene family in ^S*. 
miltiorrhiza has been previously characterized by our research 
group [31]. 

RDRs are core components of various gene silencing pathways 
and play vital roles in plant development and antiviral defense 
through regulating gene expression at the transcriptional and post- 
transcriptional levels [52]. RDR genes have been identified in 
various plants, such as Arabidopsis [14], rice [16] and maize [17]. 
However, many of them were computationally predicted only. For 
instance, among the six Arabidopsis AtRDRs, only three, including 
AtRDRl, AtRDR2 and AtRDRG, were cloned [7]. Of the five rice 
OsRDRs, only OsRDR6 were experimentally isolated [53]. 
Through genome-wide analysis of the working draft of the S. 
miltiorrhiza genome, we identified a total of five SmRDRs in S. 
miltiorrhiza. The full-length cDNAs of all predicted SmRDRs were 
then cloned and characterized. To our best knowledge, it is the 
first set of full-length SmRDR cDNAs from S. miltiorrhiza. The 
results provide useful information for further elucidation of RDR 
functions in S. miltiorrhiza. 

Using a comprehensive approach, which combines sequence 
feature, gene structure, conserved domain and phylogenetic 
analyses, we characterized the five identified SmRDRs. The 
genomic sequence of SmRDRS generating the longest SmRDR 
cDNA contains three introns, all of which locate in the open 
reading frame (ORF); whereas the SmRDR5 gene that produces 
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Figure 8. 5a77/?Z?/?5 responsive to CMV infection. The expression patterns were analyzed using the quantitative RT-PCR method. PCR was carried 
out in triplicates for each biological sample. Three independent biological replicates were performed. SmUBQIO was used as a reference. Fold 
changes of SmRDRs in leaves of 5. miltiorrhiza plantlets infected with CMV for 12, 24, 48 and 72 h are shown. The level of transcripts in leaves 
inoculated with phosphate buffered saline (CK) was arbitrarily set to 1 and the level in tissues inoculated with CMV was given relative to this. Error 
bars represent standard deviations of mean value from three biological replicates. 
doi:1 0.1 371 /journal. pone.0095117.g008 



the shortest SmRDR cDNA has seventeen introns with sixteen in 
ORF, suggesting the wide range of intron number in SmRDRs. The 
results are consistent with those from Arabidopsis and rice. The 
numbers of introns in AtRDR and OsRDR ORFs vary from 1 
[AtRDRG) to 17 {AtRDRS, AtRDR4 ^nd AtRDR5) and 1 {OsRDR6) to 
18 {OsRDR4), respectively [16]. Examination of intron numbers in 
the ORF of S. miltiorrhiza, Arabidopsis and rice RDR genes suggest 
that RDRs in each phylogenetic group have similar number of 
introns (Fig. 1) [16]. For instance, all RDR2 group members, 
including SmRDR4, AtRDR2 and OsRDR2, have 3 introns. The 
members in the RDRS group contain 16-18 introns. It suggests 
the close evolutionary relationship of RDRs in a phylogenetic 
group. 

Phylogenetic analysis of RDR proteins in S. miltiorrhiza, 
Arabidopsis and rice showed that RDRs clustered into four distinct 
groups (Fig. 5). Members of the RDRl, RDR2 and RDR6 groups 
are RDRcx-type proteins characterized by the DLDGD signature 
(Fig. 4). There is only one RDRoc-type RDR in the most recent 
common ancestor of plants, animals and fungi [47] . However, 
there are four in ^S". miltiorrhiza, of which, two are in the RDR 1 
group, one belongs in the RDR2 group, and the other one is 
included in the RDR4 group (Fig. 5). It suggests the occurrence of 
gene duplication for the RDRoc-type gene during S. miltiorrhiza 
evolution. Since these RDRcx-type RDRs clustered into different 
groups, they could be functionally diversified. All of the RDRy- 
type RDRs, which are characterized by the DFDGD signature 



(Fig. 4), clustered in the RDRS group (Fig. 5). It includes three 
Arabidopsis RDRs (AtRDRS-AtRDR5), two rice RDRs (OsRDRS 
and OsRDR4); however, there is only one ^S". miltiorrhiza RDRs 
(SmRDR5). It suggests that gene duplication has occurred for the 
RDRy-type gene during Arabidopsis and rice evolution but not in S. 
miltiorrhiza evolution. 

Gene expression analysis showed that SmRDRs were expressed 
in all of the tissues detected and responsive to MeJA treatment and 
CMV infection. All five SmRDRs were suppressed under MeJA 
treatment (Fig. 7), whereas SmRDRl, SmRDR2 and SmRDRS were 
up-regulated after CMV inoculation (Fig. 8). The involvement of 
S. miltiorrhiza SmRDRs in anti-viral defense is consistent with 
previous results for Arabidopsis AtRDRl, an ortholog of SmRDRl 
and SmRDR2, and AtRDRG, an ortholog of SmRDRS [7]. Further 
genetic manipulation of SmRDRs will shed light on the small RNA 
biogenesis pathways and RDR functions in S. miltiorrhiza. 

Supporting Information 

Table SI Primers used for 5 -RACE of SmRDRs. 

(DOC) 

Table S2 Primers used for 3 -RACE of SmRDRs. 

(DOC) 

Table S3 Primers used for amplification of full-length 
SmRDR cDNAs. 

(DOC) 
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